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I.  INTRODUCTION 


Current  Army  investigation  of  projectiles  containing  chemical  payloads 
and  smoke/incendiary  agents  has  generated  new  interest  in  liquid-filled 
shell.  These  often  exhibit  flight  behavior  different  from  that  of  solid- 
payload  projectiles  because  of  an  overturning  moment  exerted  by  the  liquid  on 
the  casing.  In  fact,  this  liquid  moment  has  produced  instabilities  which  are 
absent  for  solid-payload  shell. 

The  work  of  Stewartson1,  which  assumed  inviscid  flow,  and  the  viscous 
corrections  of  Wedemeyer2*3  have  provided  the  basis  of  the  theoretical  predic¬ 
tions  made  heretofore.  This  theory  applies  to  liquids  in  solid-body  rotation 
in  projectiles  flying  with  small  yaw.  It  demonstrates  the  existence  of  oscil¬ 
lations  in  the  liquid  and  predicts  possible  instability  when  resonance  occurs; 
i.e.,  when  the  frequency  of  angular  motion  (nutation)  of  the  shell,  t,*  is 
equal  to  a  liquid  ei genf requency  (natural  frequency  of  free  oscillation),  Cr. 

The  Stewartson-Wedemeyer  theory  yields  values  for  several  liquid  parame¬ 
ters  in  the  region  of  resonance:  frequency  and  damping  of  free  oscillations, 
and  liquid  pitch  (in  plane)  and  yaw  (out  of  plane)  moments.  The  predictions 
of  this  theory  have  been  generally  substantiated  by  gyroscope  experiments4. 

The  Stewartson-Wedemeyer  theory  assumes  that  the  timewise  variation  of 
the  flow  variables  everywhere  in  the  container  is  the  same  as  that  of  the 
nutational  motion.  This  assumption,  applicable  for  stability  studies,  would 
only  be  valid  for  actual  shell  late  in  the  flight,  after  the  fluid  has  been 
fully  spun-up  and  coned-up;  i.e.,  liquid  transients  have  decayed.  This 
assumption  will  be  retained  in  the  present  work,  yielding  a  time-independent 
problem. 


*  Definitions  are  given  in  List  of  Symbols ,  p.  41. 

1.  K.  Stewartson ,  "On  the  Stability  of  a  Spinning  Top  Containing  Liquid 
Journal  of  Fluid  Mechanics ,  Vol.  5 ,  Part  4,  September  1959,  pp.  577-592. 

2.  E.  H.  Wedemeyer ,  "Dynamics  of  Liquid-Filled  Shell:  Theory  of  Viscous 
Corrections  to  Stewartson 's  Stability  Problem ,"  BEL  Report  1287 ,  June 
1965.  AD  472474. 

3.  E.  H.  Wedemeyer ,  "Viscous  Corrections  to  Stewartson' s  Stability  Criterion," 
BEL  Report  1325,  June  1966.  AD  489687. 

4.  R.  Whiting  and  N.  Gerber,  "Dynamics  of  a  Liquid-Filled  Gyroscope:  Update 
of  Theory  and  Experiment , "  BRL  Technical  Report  ARBRL-TR-02221,  March 
1980.  AD  A083886. 
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Previously  a  theory5  was  developed  that  uses  the  linearized  Navier-Stokes 
equations  to  obtain  the  free  oscillations  of  the  fluid.  A  modal  analysis 
resulted  in  an  eigenvalue  problem  for  Cp  and  the  decay  rate  of  the  waves, 

Cj.  The  complexity  of  this  eigenvalue  problem  requires  a  large  scale  computer 

program.  A  major  part  of  this  program  is  also  required  in  the  present  problem 
which  studies  the  forced  oscillations  resulting  from  projectile  nutation.  The 
results  should  be  more  accurate  than  those  of  the  Stewartson-Wedemeyer  theory, 
which  employs  inviscid  flow  equations  plus  viscous  corrections  at  sidewall  and 
endwal Is. 

This  work  addresses  the  task  of  calculating  liquid  moment  by  an  extension 
of  the  above-mentioned  procedure5.  The  method  is  to  solve  the  internal  flow 
problem,  calculate  the  liquid  moment  by  integrating  the  computed  pressure,  and 
then  obtain  the  resultant  motion  from  the  dynamical  equations  of  yaw.  This 
study  will  treat  only  the  case  of  the  completely-filled  shell.  The  moment  due 
to  viscous  shear  will  be  reported  separately.  The  approach  is  similar  to  that 
of  Reference  5;  namely,  to  apply  a  modal  analysis  (separation  of  variables)  in 
the  flow  solution  which  gives  rise  to  ordinary  differential  equations,  and  in 
the  process  making  a  correction  required  to  compensate  for  neglect  of  the  no¬ 
slip  conditions  at  the  endwalls  in  the  modal  analysis. 

The  results  of  the  present  analysis  will  be  compared  with  experimental 
results  whenever  possible  and  with  theoretical  results  of  Murphy6.  The  latter 
uses  Wedemeyer's  viscous  corrections  at  both  sidewall  and  endwalls  to  resolve 
Stewartson's  equations. 

Some  nomenclature  that  has  been  used  regularly  in  the  past2"5'7  will  now 
be  supplanted  by  nomenclature  from  Reference  6. 

II.  EQUATIONS  OF  YAWING  MOTION 

We  shall  limit  our  consideration  to  straight  trajectories  and  small  yaw 
angles.  We  introduce  two  coordinate  systems.  The  first  is  an  inertial, 
earth-fixed  system  of  axes  x,  y,  z.  The  x-axis  coincides  with  the  projectile 
velocity  vector,  and  the  z-axis  lies  in  the  vertical  plane;  then  the  y-axis  is 
directed  so  as  to  form  a  right-handed  system.  The  second  system  is  the  aero- 

ballistic  x,  y,  z  non-rolling  system  which  has  the  x-axis  along  the  projectile 
axis  of  symmetry  and  the  z-axis  initially  in  the  vertical  plane.  These  sys¬ 
tems  are  shown  in  Figure  1;  the  y  and  z  axes  are  omitted  for  clarity. 


5.  C.  W.  Kitchens,  Jr.,  N.  Gerber,  and  R.  Sedney,  " Oscillations  of  a  Liquid  in 
a  Rotating  Cylinder:  Part  I.  Solid-Body  Rotation,"  BRL  Technical  Report 
ARBRL-TR- 02081,  June  1978.  AD  A057759. 

6.  C.  H.  Murphy,  "Angular  Motion  of  a  Spinning  Projectile  with  a  Viscous 
Liquid  Payload,"  BRL  Technical  Report  in  preparation. 

7 .  Engineering  Design  Handbook ,  Liquid-Filled  Projectile  Design ,  AMC  Pamphlet 
706-165,  April  1969.  AD  853719. 
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The  x  =  0  and  £  =  0  values  are  located  at  the  midplanes  of  the  unyawed 
and  yawed  cylinders,  respectively.  The  5f-axis  is  nutating  about  the  x-axis  at 
the  angle  Kj  =  K^(t).  The  components  of  the  projection  in  the  y,  z  plane  of  a 

unit  vector  lying  on  the  £-axis  are  denoted  by  nY£  and  n^. 

The  yawing  motion  is  characterized  by  two  variables,  3  and  g.  The  angle 
of  attack,  a,  in  the  aerobal 1 i Stic  system  is  the  angle  in  the  vertical  plane 
measured  also  from  the  x-axis  to  the  velocity  vector;  the  angle  of  slideslip, 
3,  is  the  angle  in  the  horizontal  plane  measured  also  from  the  £-axis  to  the 
velocity  vector.  For  the  small  yaw  angles  considered,  3  s-n^  and  g  ~~nYE* 

It  is  convenient  to  combine  3  and  g  into  a  single  complex  variable: 

5  =  g*-i  3  -  -(nY£  +  i  (1) 

The  fluid  pressure  forces  on  the  cavity  surfaces  produced  by  the  motion 
give  rise  to  a  moment  on  the  projectile.  The  spin-decelerating  component, 

Ml£,  is  zero  here;  the  other  components  can  be  represented  in  complex  form, 

Mly~  +  iM|_£.  We  shall  consider  only  the  liquid  moment  acting  on  the  projec¬ 
tile.  The  resulting  model  is  adequate  for  comparing  theoretical  outputs;  the 
liquid  moment  can  be  added  to  the  other  moments  acting  on  shell  or  gyroscrope 
as  required.  The  differential  equation  of  yawing  motion  is  * 


I  d2?/dt2  -  i+Ix  d?/dt  +  I  M?  =  i  (MLf  +  1Ml?).  (2) 

The  quantity  Ix  is  the  moment  of  inertia  of  the  empty  axisymmetric  shell  about 

its  longitudinal  axis.  \y  is  the  traverse  moment  of  inertia  of  the  empty 

shell  about  its  center  of  gravity.  The  spin  rate  of  the  shell  is  <f>  (taken  to 

be  positive);  t  is  time.  The  term  IyMc  is  an  aerodynamic  moment  for  a 

projectile.  For  a  gyroscope  this  term  is  a  gravitational  moment  arising  from 
the  separation  of  center  of  gravity  and  pivot  point,  and  in  most  experiments 
is  zero. 

In  general  there  is  an  interaction  between  the  motion  of  the  projectile 
and  the  liquid  motion.  Here  we  shall  specify  the  motion  of  the  projectile. 

In  particular  the  cylinder  is  nutating  with  constant  frequency  and  exponen¬ 
tially-growing  yaw: 


c=  (K„e”#)e1(T*)  =  Kje1*1  . 


(3) 


*  This  is  Eq. 
hand  side. 


(2.4)  of  Reference  6  with  only  the  liquid  moment  on  the  right- 
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where 


=  K0eeT^* 


(4) 


4>1  =x$t ,  f = ( 1  -  i  e )  x  • 


Here  K0  is  the  magnitude  of  the  yaw  at  time  t  =  0;  x  is  the  nutational  fre- 

* 

quency  divided  by  <$>,  and  the  yaw  grows  when  ex  >  0.  Also  is  the  yaw  ampli¬ 
tude,  and  is  the  angular  orientation*  of  the  X-axis  in  the  x,  y,  z  system 
as  shown  in  Figure  1. 


The  motion  of  the  projectile  enters  the  flow  problem  via  the  boundary 
conditions.  Under  the  assumption  that  the  flow  is  in  phase  with  the  motion  of 
the  shell,  the  pressure  disturbance  will  have  the  time  dependence  of  Eq.  (3), 
and  consequently  the  liquid  moment  will  also  have  this  form.  A  nondimensional 
liquid  moment  coefficient,  C^,  is  now  defined**: 


Mj-Y  +  i  =  x  C|_f^  K^o 


Hi 


(5) 


where  mL  is  the  mass  of  the  liquid.  is  a  complex  quantity  whose  real  part 

represents  a  moment  that  changes  the  yaw  angle,  and  whose  imaginary  part 
changes  the  nutation  rate.  Thus: 


CLM  E  CLSM  +  1  CLIM»  (6) 

where  CLsm  and  CLjM  represent  the  "liquid  side  moment"  and  "liquid  in-plane 
moment",  respectively.  That  Cj_5M  and  represent  these  moments  can  be 

demonstrated  by  taking  the  scalar  products  of  the  moment  vector  (MLy,  M^) 
with  a  unit  vector  parallel  to  the  yaw  vector,  (cos  ,  sin  and  a  unit 

vector  normal  to  the  yaw  vector,  (-  sin  cos  ,  respectively. 

III.  FLOW  PROBLEM 

A.  Flow  Equations  and  Boundary  Conditions 

Cylindrical  polar  coordinates  are  introduced  in  the  earth-fixed  frame 
y  =  r  cos  9,  z  =  r  sin  9,  x  =  x,  (7) 


*  For  simplicity  the  angle  of  attack  is  assumed  to  he  initially  zero  and  the 
angle  of  sideslip  to  be  initially  positive ;  i.e.>  <J>iq  of  Reference  6  is  zero. 

** See  Eq.  (2.7)  in  Reference  6. 
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and  in  the  non-rotating  aeroball istic  frame 


y  =  r  cos  e,  z  =  r  sin  8,  x  =  x.  (8) 

All  the  above  terms  are  non-dimensional;  lengths  and  distances  are  non- 
dimensional  ized  by  a,  the  cross-sectional  radius  of  the  cylinder.  The  rela¬ 
tionship  between  the  two  sets  of  cylindrical  coordinates  for  small  K0, 

obtained  from  the  orthogonal  transformation  between  x,  y,  z  and  x,  y,  z,  is 
r  =  ?  -  Kj  (?-£)  cos  Uj-0)  +  0  (K02)  (9a) 


8=9-^  [(x- t)/r]  sin  (<j>1-e)  +  0  (K02) 


(9b) 


x  =  x  +  K1  r  cos  (<^-6)  +  0  (K02), 


(9c) 


where  i  is  the  non-dimensional  x-coordinate  of  the  pivot  point. 

To  obtain  the  flow,  we  assume  a  small  disturbance  to  a  basic  flow,  which 
is  taken  to  be  solid-body  rotation  in  an  unyawed  cylinder.  The  Navier-Stokes 
equations  are  linearized  to  produce  the  perturbation  equations.*  The  flow 
variables  are  the  radial,  azimuthal,  and  axial  velocity  components,  and  pres¬ 
sure,  given  here  in  non-dimensional  form: 

u  =  U  -  K0u ,  v  =  V  -  K0v,  w  =  W  -  Kqw,  p  =  P  -  K0p.  (10) 

The  symbols  u,  v,  w,  and  p  represent  the  total  values.  U,  V,  W,  and  P  are  the 

basic  undisturbed  variables;  and  u,  v,  w,  and  p  are  perturbation  variables** 
of  order  one.  For  solid  body  rotation  the  basic  flow  is 

U  =  0,  V  =  r,  W  =  0,  3P/8r  =  r.  (11) 

The  velocity  components  are  non-dimensional  ized  by  a$,  and  pressure  by  pa  24>z , 
where  p  is  the  density  of  the  liquid. 

On  substituting  Eqs.  (10)  and  (11)  into  the  Navier-Stokes  equations  and 
linearizing  with  respect  to  K0,  we  obtain  the  perturbation  equations  in  non- 

dimensional  form: 


*  These  ave  Eq.  (3)  in  Reference  5. 

**  The  negative  signs  in  Eq.  (10)  were  employed  to  comply  with  the  nomencla¬ 
ture  of  Reference  6. 
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★  ★ 
ve  +  rwx  = 


+  u0  -  2v  =  -pr  +  Re'1  [v2u  -  r-2  u  -  2r'2  v0] 

+  Vq  +  2u  =  -r'1  p0  +  Re-1  [V2v  -  r-2  v  +  2r"2  u0] 
w£  +  w0  =  -px  +  Re-1  v2w, 

where  subscripts  denote  partial  differentiation,  t  =  $t,  and 
v2  e  a2/ ar2  +  r_1  a/ ar  +  r-2  a2/ae2  +  a2/ax2. 

The  Reynolds  number  is 
Re  =  $  a2/v, 

where  v  is  the  kinematic  viscosity  of  the  fluid. 


(12a) 
(12b) 
(12c) 
( 1 2d ) 

(13) 

(14) 


The  boundary  conditions  are:  no  flow  through  the  bounding  walls  and  no 
slip  along  them;  i ,e. , 

Q(r*=l)  =  »(P=1)  =  0,  9(r*=l)  =  1  (side)  (15a) 

0(x=±c)  =  w(x=±c)  =  0,  v(x=±c)  =  ?  (end)  (15b) 

where  u,  9,  w  are  non-dimensional  cylindrical  velocity  components  in  the 

••v  S 

aeroballistic  system  r*,  e,  x;  and  c  is  the  half-height,  c,  of  the  cylinder 
divided  by  a. 


The  boundary  conditions  must  be  transformed  to  the  coordinates  used  in 
Eq.  (12);  the  transformations  of  Eq.  (9)  are  used  to  accomplish  this.*  The 
resulting  non-homogeneous  boundary  conditions  in  the  inertial  system  are: 


“wall  =  Rea1  Mn-O  exP  ( i  (f $t-e ) } ]  +  0(Ko)  (16a) 

vwail  =  -(x-i)  Real  [(!-f)  exP  {i(f$t-e)}]  +  0(Ko)  (16b) 

wwan  =  -r  Real  [-i(l-f)  exp  (i(f$t-e)}]  +  0(Ko).  (16c) 


(Recall  that  f  =  (l-ie)x.) 


*  A  derivation  is  outlined  in  Appendix  <4. 
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The  t  and  e  dependence  of  the  boundary  conditions  is  satisfied  by  a  flow 


solution  of  the  form 

u  =  Real  (uq)  =  Real  [u(r,x)  exp  (i ( f $t-0 ) } ]  (17a) 
v  =  Real  (*q)  =  Real  [v(r,x)  exp  (i(f$t-6)}]  (17b) 
w  =  Real  (Sc)  =  Real  [w(rsx)  exp  (i(f$t-e)}]  (17c) 
P  =  Real  (pc)  =  Real  [p(r,x)  exp  { i ( f^t-9 ) } ] »  (17d) 


where  u,  v,  w,  and  p  are  complex  functions.  The  functions  u^,  *q,  w c,  and  pq 

are  clearly  also  solutions  of  Eq.  (12).  We  substitute  these  complex  solutions 
into  Eq.  (12)  to  obtain: 

ru  +  u  -  iv+  rw  =  0  (18a) 

i ( f - 1 ) u  -  2v  =  -p  +  (1/Re)[u  +  u  /r  -  2u/r2 

-r  -rr  -r 

+  u  +  2iv/r2]  (18b) 

-xx 

i ( f-1 ) v  +  2u  =  ip/r  +  (1/Re)[v  +  v  /r  -  2v/r2 

-  -rr  -r 

+  v  -  2iu/r2]  (18c) 

-xx 

i(f-l)w  =  -p  +  (1/Re)[w  +  w  /r  -  w/r2  +  w  ].  (18d) 

-x  -rr  -r  -  -xx 

The  wall  boundary  conditions  of  Eq.  (16)  become 

“wall  ■  (X-‘>'  V-wall  =  -(1-f)  (X-‘)’  (19) 

"wall  ■  U1-f>  r ■ 

Axial  boundary  conditions  are 

u(r=0)  -  iv(r*0)  =  0,  w(r=0)  =  0,  (20) 

p(r=0)  =  0. 
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These  are  required  by  the  equations  themselves,  assuming  that  all  variables 
and  their  derivatives  are  finite  at  r=0;  a  derivation  is  given  in  Appendix  B. 


The  solution  to  Eq.  (18)  is  expressed  as  the  sum  of  two  solutions: 


u=u+u,  v  =  v  +  v  (21) 

-p  -H  -  -p  -H 

w  =  w  +  w  ,  p  =  p+p, 

-p  -H  -  -p  -H 

where  the  particular  solution  (u  ,  v  ,  w  ,  p  )  is 

-p  -p  -p  -p 

u  =  -i[(l-f)2/(l+f)]x  +  i  (l-f)i.  (22a) 

-P 

v  =  -[(l-f)2/(l+f )]x  +  (l-f ) A  (22b) 

-P 

w  =  i (l-f)r  (22c) 

~P 

p  =  -(l-f)2rx+(l-f2)  ir.  (22d) 

-P 

Sidewall  boundary  conditions  are  now: 

u  (r=l)  =  -i[2f (1-f )/ (l+f)]x  (23a) 

-H 

v  (r=l)  =  -[2f (1-f )/(l+f)]x  (23b) 

-H 

w  (r=l)  =  0.  (23c) 

-H 

Endwall  boundary  conditions  are: 

u  (x=c)  =  -i [2f (1-f )/(l+f) ]c ,  u  (x=-c)  =  -u  (x=c)  (24a) 

-H  -H  -H 

v  (x=c)  =  -[2f (1-f )/ (1+f )]c,  v  (x=-c)  =  -v  (x=c)  (24b) 

-H  -H  -H 

w  (x=c)  =  w  (x=-c)  =  0.  (24c) 

-H  -H 

Axial  boundary  conditions  are  the  same  as  those  of  Eq.  (20): 

-H(r=°H-H(r=0)  -  wH(r=0)  =  pH(r=0)  =  0.  (25) 


The  total  problem  is  now  finally  given  by  Eqs.  (18),  (23),  (24),  and  (25) 

for  u  ,  v  ,  w  ,  and  p  . 

-H  -H  -H  -H 
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B.  Modal  Analysis:  Separated-Variable  Solutions 
1.  Form  of  Solution. 


The  solutions  to  Eq.  (18)  will  be  obtained  using  separation  of  vari¬ 
ables  with  the  corrected  endwall  condition  described  in  the  next  section. 
This  gives  rise  to  a  non-standard  eigenvalue  problem  for  the  x  variation. 
Separation  of  variables  requires  solutions  of  the  form 


uH(r,x)  =  Rj(r)  X^x)  w^r.x)  =  Mr)  X3(x)  (26) 

vH(r,x)  =  R2(r)  Mx)  PH(r»x)  =  R4^  *4^ 

Substituting  these  into  Eq.  (18)  yields  ordinary  differential  equations  for 
the  R's  and  X's.  It  can  be  shown,  see  Appendix  C,  that  the  X.-  must  satisfy 

the  harmonic  equation 


d  2X.j/dx  2  +  A2Xi  =  0  (27) 


for  i  =  1,  2,  3,  4,  where  x  is  the  eigenvalue.  The  equations  for  R-j  are  given 
in  Eq.  (33)  below.  The  conclusion  is  that  the  general  form  (26)  reduces  to 
that  assumed  in  the  modal  analysis  of  Reference  5. 

Another  possible  approach  is  the  method  of  Hall8,  which  was  developed  for 
circular  Couette  flow  in  a  finite  length  cylindrical  annulus.  It  would  use  a 
modified  form  of  the  theory  of  Reference  5,  and  the  no-slip  condition  on  the 
endwall  would  be  satisfied  at  a  finite  number  of  points.  One  feature  of  this 
method,  in  common  with  the  approach  adopted  here,  is  that  the  x^  (or  k  in 
Reference  5)  are  complex. 

The  separation  of  variables  leads  to  the  result  that  all  flow  variables 
have  the  form 


R(r)[Aj  sin  Xx  +  A2  cos  Xx], 


The  sidewall  boundary  conditions  Eqs.  (23a)  and  (23b),  show  that  u  and  v  are 

-H  -H 

odd  functions  of  x  at  r=l;  thus  series  expansions  of  eigenfunctions  should 
contain  only  the  odd  functions;  i.e.,  the  sin  xx's  for  u^  and  v^.  Inspection 

of  Eqs.  (18a)  and  (18b)  containing  these  u  and  v  solutions  shows  that  w 

-H  -H  -H 


8.  P.  J.  Blennerhasset  and  P.  Hall ,  "Centrifugal  Instabilities  of 

Circumferential  Flow  in  Finite  Cylinders:  Linear  Theory Proa.  Roy.  Soa. 
London  A-365,  pp.  191-207,  1979. 
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must  then  contain  only  cos  ax  and  p  must  contain  only  sin  Ax.  Thus  we  assume 

-H 

a  solution  in  the  form  of  the  following  infinite  series:* 


-H  =  k?i  Uk(r)  $in  XkX’ 

v  =  S  v(r)  sin  Al,x 
-H  k=l  k 

(28) 

-H  =  'k?i  Wk(r)  C°S  XkX’ 

-H  =  Pk(r)  Sin  XkX* 

A  A  A  A 

The  u^,  v^,  w^,  are  complex  functions  which  are  solutions  of  Eq.  (33) 

below.  In  the  procedure  described  next,  for  satisfying  approximately  the  end- 
wall  boundary  conditions,  Eq.  (24),  the  eigenvalues  x^  are  complex;  the  "un¬ 
corrected"  form  of  this  boundary  condition  would  give  real  and,  therefore, 
the  normal  modes. 

2.  Corrected  Endwall  Boundary  Condition 

The  above  solution  cannot  satisfy  all  the  endwall  boundary  condi¬ 
tions,  Eq.  (24).  For  xk  =  k7r/(2c),  where  k  is  odd,  it  satisfies  Eq.  (24c)  but 

fails  to  satisfy  the  no-slip  condition.  To  correct  for  this,  a  procedure 
analogous  to  that  of  Wedemeyer2*3  is  adopted.  The  technique  of  matched 
asymptotic  expansions  is  used  to  derive  the  correction;  it  is  explained  in 
detail  in  Appendix  D.  The  three  boundary  conditions  of  Eq.  (24)  are  required 
to  determine  the  solution  to  Eq.  (18).  It  is  shown  in  Appendix  D  that  the 
outer  solution  to  Eq.  (18)  is  determined  by  only  one  boundary  condition  on 
w^.  The  outer  solution  is  matched  to  the  inner,  boundary  layer  type,  solution 

which  satisfies  the  no-slip  condition.  The  one  condition  on  the  outer  solu¬ 
tion  for  w^  can  have  different  forms  for  various  degrees  of  approximation. 

Thus, 


wH  =  0  at  x  =  t  c 


gives  the 

solution  to  0  (Re-1/2); 

this  is  the  uncorrected  solution  with 

real  xk. 

The  next  approximation 

w^  -  6c  ^h/3x  =  0 

at  x  =  c 

(29a) 

w^j  +  6c  ^h/3x  =  0 

at  x  =  -c 

(29b) 

gives  the  solution  to  0  (Re-*);  this  is  the  corrected  solution  with  complex 

V 


*  The  negative  sign  appears  in  the  w  series  for  the  purpose  of  facilitating 

~H 


the  adaptation  of  existing  programming  to  the  present  problem. 
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Here  non-dimensional  6c  is  given  by  the  following  sequence: 


a  =  2 _1/2  Re1/2(l-i)(3-f)1/2  (30) 

3  =  2~1/2  Re1/2(l+i)(l+f)1/2. 


The  complex  square  roots  are  chosen  to  be  the  ones  that  make  the  real  parts 
of  a  and  3  positive.  Detailed  expressions  are  found  in  Eq.  (A. 5.)  of  Refer¬ 
ence  5. 

The  theory  giving  conditions  (29)  is  asymptotic  for  Re  +  ®.  The  lower 
limit  on  Re  that  will  give  a  certain  accuracy  in  the  solution  must  be  deter¬ 
mined  by  comparison  of  the  results  with  experiment  or  by  getting  the  solution 

using  the  next  term  for  the  boundary  conditions  to  0  (Re-3^2). 

Each  cos  xkz  term  of  w  in  Eq.  (28)  must  individually  satisfy  Eq.  (29). 

=  -H 

Substitution  at  x  =  ±  c  yields  the  functional  equation  for  the  denumerable  set 
of  permissible  Ak's, 


cos  Akc  +  Xk  6c  sin  xkc  =  0.  (31) 

For  | 6c | /c<<l , 

Ak  s  K=  (kir)/[2(c-6c)]  (k  odd)  (32) 

3.  Ordinary  Differential  Equations  for  Radial  Variation 

The  solution  in  Eq.  (28)  portrays  the  modal  composition  of  the 
perturbed  flow  in  the  axial  direction;  each  term  of  the  four  series  is  a 
solution  to  Eq.  (18).  When  the  four  k  th  terms  are  substituted  into  Eq.  (18), 
the  sin  Akx  and  cos  Xkx  terms  cancel  out,  leaving  the  following  set  of  ordi¬ 
nary  differential  equations  for  uk,  vk,  wk,  and  pk  (where  '  =  d/dr),  omitting 
subscripts  k  on  the  variables: 


ru  ’  +  u  -  iv  +  Akrw  =  0 


Re-1  u"  +  (Re  r)"^1  +  [i  (1-f)- Re ~ 1  ( 2 / r 2  +  Xk2)]u  + 


[2+2i (Re  r2)-1]v  =  p  ' 


(33a) 

(33b) 
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(33c) 


Re-1  v  "  +  (Re  r)~*  v  '  +  [i  (l-f)-Re-1(2/r2+Ak2)]v  - 

A  At 

[2+2i(Re  r2)_1]u  =  -i  p/r 

A  A  A  At 

Re-1  w"+  (Re  r)-1  w '+[i  (l-f)-Re_1(l/r2+Xk2)]w  =  -Xkp. 


(33d) 


These  equations  are  converted  to  canonical  form  in  order  to  be  integrated 
numerically;  i.e., 


yi  '  =  dyi/dr  =  f^r,  yls  y2,...y6), 


where 


Yl  =  u 

A  At 

y2  =  u  -  iv 


y3  =  v 

A 

Y4  =  w 


i  =  1,  2,. ..6, 

y5  =  w' 

A 

y6  =  p* 


(34) 


After  the  required  manipulations  are  performed,  the  following  sixth  order 
system  is  obtained: 


Yl'  =  -(y2/r)  -  Ak  y4 
y2'  =  -(y2/r )  -  1  Y3“xk  Y4 


(35a) 

(35b) 


y3 '  =  2(Re  +  i  r"2^  +  i  (P+r*2)  (y2-y3)  -  y3/r  -  i(Re/r)y6  (35c) 


where 


y4  =  Y5 

y5'  =  e  y4  -  y5/r  -  xk  Re  y6 
y6'  =  -Re'lfl  yx  +  [2i-(Re  r2)"l](y2-y1)  + 
i (Re  r)-1y3-Ak  Re_1y5 

3  =  r-2  +  Xk2  -i  Re  (1-f). 


(35d) 
(35e) 
( 35  f ) 


(36) 


4.  Boundary  Conditions  for  Radial  Equations 

There  are  three  boundary  conditions  at  r=0  and  three  at  r=l.  As  a 
consequence  of  Eqs.  (25)  and  (28),  conditions  at  r=0  are  (reintroducing  the 


A  A 


index  k)  uk-ivk  =  wk  =  pk  =  0,  or 

Y2k(0)  =  Y4k(0)  =  Y6k(®)  =  0* 

From  Eqs.  (23)  and  (28)  we  obtain  at  r  =  1: 


(37) 
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(38) 


-i[2f(l-f)/(l+f)]  x  =S  uk(l)  sin  Akx 

In 

-  [ 2f (l-f)/(l+f)]  x  vk(l)  sin  *kx 

y  * 

0  =  -  y  wk(l)  cos  Akx. 

A  A 

In  order  to  obtain  ,uk(l)  and  vk(l),  we  expand  the  function  x  in  the  inter¬ 
val  -c<x<c  in  a  series  in  Xk(x)  =  sin  xkx: 


sin  Akx. 


(39) 


The  functions  Xk(x)  are  the  eigenfunctions  of  the  x-eigenvalue  problem  dis¬ 
cussed  in  Section  III.B.l.  When  6c  =  0  the  xk  are  real  and  the  eigenvalue 
problem  is  a  standard  self-adjoint  problem.  The  eigenfunctions  are  orthogonal 
and  form  a  complete  set;  in  fact,  the  bk  are  the  Fourier  coefficients.  For 

6c*0,  and  since  6c  is  complex,  a  standard  self-adjoint  problem  is  not  ob¬ 
tained.  The  eigenfunctions  are  not  orthogonal.  They  are  biorthogonal  with 
respect  to  the  solutions  of  the  adjoint  problem,  which  property  enables  the  bk 

to  be  determined.  Some  of  the  details  are  discussed  in  Appendix  E;  here  we 
merely  state  the  coefficients: 


(2/Ak2)  [1  +  (Xk6c)2]  sin  Akc 
C  [1  +  (Ak  6c)2]  -  6c 


(40) 


From  Eqs.  (34),  (38),  and  (39)  the  sidewall  boundary  conditions  are: 


yikd) 

=  uk(l)  =  -i  [2f(l-f)/(l+f)]  bk 

(41a) 

Y2k(1) 

A  A 

=  uk(D  "  1  vk(!)  =  0 

(41b) 

y4k ( 1 ) 

• 

o 

II 

r- H 

<  5 

Ml 

(41c) 

An  alternative  set  of  coefficients,  ak,  are  also  used  for  the  sin  xRx  expan¬ 
sion  of  x;  actually,  these  had  been  used  before  the  bk  formula  of  Eq.  (40)  was 
derived.  They  are  determined  by  optimizing,  in  a  least  squares  sense,  the  fit 
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kf 

of  the  partial  sum,  a^  sin  x  (k  odd),  to  the  function  x  over  the  inter- 

k=l 

val  -c<x<c.  Thus,  the  (kf-l)/2  a^'s  are  calculated  which  minimize  the  follow¬ 
ing  integral  (where  k  is  odd): 

fc  kf 

I  (a^  a3,...akf)  =  1  _  I  x  -  £  ak  sin  Xkx  I2  dx*  (42) 

J-c  k=l 

In  our  work  the  minimization  was  carried  out  with  a  preprogrammed  computer 
routine  that  involved  an  iterative  operation.  The  Fourier  series. 


x  = 


L)/2  8c/ (nk)2 


sin 


kir/ (2c)  x 


) 


(k  odd), 


furnished  good  initial  guesses  for  the  a^'s.  We  have  generally  used  eight 
terms  in  the  series  of  Eq.  (39)  or  its  counterpart  with  a^'s.  These  have 
furnished  sufficiently  accurate  representations  of  the  function  x. 


5.  Numerical  Procedures. 


The  differential  system,  Eq.  (35),  and  initial  conditions,  Eq.  (37), 
constitute  the  same  numerical  integration  problem,  including  orthonormaliza¬ 
tion,  as  does  the  eigenf requency  determination5.  We  refer  the  reader  to 
Reference  5  for  a  detailed  description  of  the  numerical  procedures.  One 
difference  arises  because  the  terminal  conditions,  Eq.  (41),  are  non-homogene- 
ous. 


In  the  orthonormalization  process5,  the  interval  Ocrcl  is  divided  into  N 
sub-intervals.  In  each  subinterval  three  linearly  independent  solutions  for 

u,  v,  w,  p  which  originally  satisfy  the  initial  conditions  at  r  =  0  are  found, 

and  these  are  combined  linearly  so  that  the  total  solution  is  continuous  at 
the  boundaries  of  the  sub-intervals.  The  process  of  determining  the  coeffi¬ 
cients  of  the  combinations  begins  in  the  last  (N  th)  subinterval  by  applying 
the  inhomogeneous  sidewall  conditions  at  r  =  1  for  a  given  k: 


C1N  lJ  1 N ( 1 )  +  C2N  U2N^)  +  C3N  U3N^) 
C1N  v1n(X)  +  C2N  V2N^ 1  ^  +  C3N  V3N^) 
C1N  W1N(1)  +  C2N  W2N^ 1 )  +  C3N  w3n(:) 


u(l) 

v(l) 


0. 


(43) 
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The  quantities  on  the  right-hand  sides  are  obtained  from  Eq.  (41).  The  other 
sets  of  c^,  C2»  C3  for  the  succeeding  subintervals  can  now  be  determined  in 

sequence  proceeding  from  r  =  1  to  r  =  0  by  the  method  described  in  Section 
III.E.  of  Reference  5. 


IV.  EVALUATION  OF  PRESSURE 

The  flow  variable  of  primary  interest  here  is  the  pressure.  From  Eq. 
( 17d) , 

p  =  ee*Tt  [-pi  sin  (x$t  -  9)  +  pR  cos  (x*t  -  9)J  +  0  (K0) ,  (44) 


where 


P  (r,x)  =  pR  +  i  p  .  (45) 

The  total  pressure,  to  within  an  additive  constant,  recalling  Eqs.  (10)  and 
(11),  is 

P  =  \  rZ  "  *0  P*  (46) 

The  liquid  moment  calculations  call  for  pressure  evaluations  at  constant 

values  of  r  and  x.  Accordingly,  the  —  r2  term  in  Eq.  (46)  is  replaced,  from 

2 

Eq.  (9a),  by 

.1-  r2  =  -L  r2  -  K0  ee*xt  r(x-fc)  cos  (x<frt  -  0)  +  0  (K02).  (47) 

Thus,  by  Eqs.  (46)  and  (47),  the  disturbance  pressure,  which  is  measured  at 
points  fixed  on  the  cylinder  surface,  is 

Dp  =  p  -  i.  r2  =  J-  r2  -  .1  r2  -  Knp. 

2  2  2  0 

Dp  =  p  -Ir2  =  -K0  eeT^  [-p  sin  (x*t  -  0)  + 

+  {p  +  r  (x-*)}  cos  (x<f>t  -  0)].  (48) 

“  R 
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The  r,  9,  x  in  the  bracketed  term  of  Eq.  (48)  can  be  replaced  by  r,  9,  x 
without  changing  the  order  of  the  approximation  since  r  =  r  +  0  (K0),  6=0+ 
0  (K0),  x  =  x  +  0  (K0).  Hence,  Dp  (r  =  1)  =  Dp  (r=l)  +  0  (K02)  and 
Dp  (x  =  t  c)  =  Dp  (x  =  ±  c)  +  0  (K02).  By  Eqs.  (21),  (22d),  and  (28), 

[? 


p  =  Real 
-R 


pk(r)  sin  X^x  +  (1-f2)  £  r  -  (1-f)2  r 


p^  =  Imag 


P|<(r)  sin  A^x  +  (1-f2)  £  r  -  (1-f)2  r  x 


The  disturbance  pressure  of  Eq.  (48)  is  rewritten  as 


=  -K0e£T^t  sin  (4>1  -  9)  +  P2  cos  (^  -  0)J 


where 


Dp 

pi  =  -  ?i- 


P2  =  p  +  r  (x-t), 
-K 


Amplitude  and  phase  of  Dp  are  denoted,  respectively,  by 


(49) 


(50) 


Cp  =  (pl2  +  P22)172  and  t  =  tan'1  (P2/pl)»  (51) 

so  that  Dp  =  -K0eex^)t  Cp  sin  j^Uj-9)  +  i|>  j . 

The  pressure  measurements  of  Whiting9  can  be  compared  with  our  pressure 
calculations.  In  these  experiments  the  liquid  is  first  completely  spun  up  in 
an  unyawed  cylinder;  then  the  cylinder  is  nutated  at  a  fixed  small  angle  of 
yaw  about  its  center  (£  =  0)  with  a  fixed  frequency.  Built-in  pressure-mea¬ 
suring  apparatus  enables  the  disturbance  pressure  to  be  measured  at  two  points 
on  the  endwall.  For  comparison  we  choose  data  from  Figures  10a,  10c,  lOd,  lOe 
in  Reference  9,  which  show  Cp  measurements  plotted  against  the  forced  coning 

frequency,  x. 

The  comparison  is  shown  in  Figures  2  and  3.  The  solid  curves  are  plots 
of  the  present  computations  of  Cp.  The  locations  of  the  peaks,  in  the  experi¬ 
mental  data  and  theoretical  results,  agree  to  within  about  2.5%.  The  largest 
discrepancies  in  amplitude  in  Figures  2a  and  3a  occur  at  the  peaks,  where  the 
differences  are  approximately  15%.  The  overall  percentagewise  agreement  is 
poorest  in  Figure  2b,  where  the  Reynolds  number  is  two  orders  of  magnitude 


9.  R.  D.  Whiting,  "An  Experimental  Study  of  Forced  Asymmetric  Oscillations  in 
a  Rotating  Liquid-Filled  Cylinder,"  BRL  Technical  Report  ARBRL-TR-02376, 
October  1981.  AD  A1 07948. 
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smaller  than  those  of  the  other  cases.  It  was  noted  above  that  the  error  in 
the  theory  increases  as  Re  decreases;  this  may  be  reflected  in  Figure  2b.  We 
note  that  the  percentagewise  scatter  in  the  measurements  is  also  greatest  for 
this  case. 

Computations  by  the  theory  of  Reference  6  are  also  plotted.  In  Figures 
2a,  3a,  and  3b  they  are  indistinguishable  from  the  solid  line  curves.  In  the 
low  Reynolds  number  case,  however,  the  differences  are  significant  and  the 
results  of  Reference  6  show  better  agreement  with  measurement.  The  eigen- 
frequency,  Cr,  computed  by  the  method  of  Reference  5,  is  indicated  for  each 

case  in  Figures  2  and  3  and  in  all  subsequent  figures.  In  diagrams  having 
peaks,  namely  those  for  pressure  coefficient,  side  moment  coefficient,  and  yaw 
growth  rate,  it  is  seen  that  the  T-value  at  the  peak  is  very  close  to  Cr  for 

large  Reynolds  numbers,  but  that  the  difference  between  these  two  frequencies 
increases  as  Reynolds  number  decreases. 


V.  LIQUID  PRESSURE  MOMENTS 


Our  objective  is  to  determine  the  moment  produced  by  the  liquid  on  the 

spinning  and  nutating  shell,  namely,  Mj_ y  +  i  Mq  of  Eq.  (5).  We  shall  evalu- 

ate  the  moment  about  the  center  of  gravity  of  the  projectile  in  the  x,  y,  z 
system.  Details  need  be  shown  for  only  one  component,  say  M^,  because  of 

axisymmetry  in  the  transverse  motion. 

MLZ  =  MLZS+  mlzt+  MLZB 

The  three  terms  on  the  right-hand  side  denote  the  moments  on  the  side,  top, 
and  bottom  walls,  respectively.* 


\1S 


(pa5*2) 


2  TT  ~  ~  ~  ~ 

p  (r=l)  (x-£)  cos  0  d9  dx 

0 


(53a) 


mLZT  =  -(pa5*2) 


(x=c)  r2  cos  9  d0 


(53b) 


^LZB  “  (pa5*2) 


.  rs.  _  ^  ~  ~  <->• 

( x=-c)  r2  cos  9  d9  dr 


(53c) 


and  are  found  in  Section  3-3.1  of  Eeference  7. 

Eqs.  (3-31)  and  (3-30a)  contain  errors  in  sign ;  there  should  he  negative  signs 
before  the  integrals. 


* Formulas  for  ^£i^s» 
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The  non-dimensional  cylinder  radius  occurring  in  the  integrand  of  Eq.  (53a)  is 
equal  to  1. 

When  p,  of  Eq.  (48),  is  substituted  into  Eq.  (53),  the  1/2  r2  term  makes 
no  contribution  to  the  integrals,  and  the  lowest  order  terms  therefore  are 

0(ko)  .  Hence  for  first  order  accuracy  in  K0  the  integrands  may  be  evaluated 


in  the  earth-fixed 

system. 

The  integrals  become: 

Mizs/(pa5$2) 

=  -IT  K1 

(sin  4>x)  fC=  -(x-t)  p  (r=l)  dx  + 

L  =  '  -c 

(54a) 

(cos  ^)  JC=  (x-t)  |pR  (r=l)  +  (x-t)}  dxj 

MLZT^pa5^2^ 

=  IT  Kj 

(sin  *!)  -r2  ^  (x=c)  dr  + 

(54b) 

(cos  ^) 

io  f2  {PR  +  r  (C"A)}  dr 

MLzB/(Pa5$2) 

=  -n  K1 

(sin  4>i)|q  -r2  Pj  (x=-c)  dr  + 

(54c) 

(cos  h)| 
«• 

*1  _  _ 

lo  r2  {?r  (x="c)  "r  (c+£) }  dr  * 

The  net  moment  on  the  endwalls  is 
MLZE  =  MLZT  +  mlzb- 


ML^E/(pa5$2)  =  Kl  j^sin  <t»i  j"0  -p2  { Pj  (x=-c)  ^  dr  + 

C0S  *1  5o  rM-R  '-R  (X="^}  +  2  rC}  dP  * 

We  observe  that  is  independent  of  i. 


(55) 
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(56) 


The  moments  of  Eqs.  (54a)  and  (55)  have  the  form: 

M|_2s=  ~  Ki  (Mjg  sin  $i  +  M2^  cos  $j)  (pa3$2) 

M|^2£=  -  Kj  (Mje  sin  4*1  +  M2j:  ^l)  (pa3^). 

The  moment  coefficients  M^,  M2$,  ^1E»  anc*  ^2E  are  functions  °f  Re»  c/a,  T» 
and  e. 


(jc= 


M1S  =  -  w  Imag|\  =  (x-£)  p  (r=l)  dx 


(57a) 


fc 


M2S  =  ir  Real  V  =  (x-t)  p  (r=l)  dx  +  (2/3)  (c3  +  3  ci2) 


(57b) 


11E 


tt  Imag  [1:  r2|p  (x=c )  -p  (x=-c  )}  dr  +  c/2  j 


(57c) 


M2e  =  ~  it  Real 


The  total  moment  is: 


ft  p2{p  (x=c)  -P  (x=-c)}dr  +  c/2  J . 


(57d) 


Mlz  =  M|_z$+  ^lz£=  (pa3$2)  (Mj  sin  $1  +  M2  cos  $j)» 


(58) 


where 


Mi  =  Mi^  +  M 


’IE* 


m2  =  m2^  +  m2e* 


(59) 


Our  computational  results  will  be  exhibited  in  terms  of  the  moment  coef¬ 
ficient,  C|_m«  defined  in  Eqs.  (5)  and  (6).  By  Eqs.  (5),  (6),  and  (58), 


Clsm  =  -M2  /  [2irrc/a],  ^^1^  =  -Mj  /  [2irrc/a], 


(60) 


Computations  of  ClsM  and  ^lim  were  performed  for  four  Reynolds  numbers  cover¬ 
ing  several  orders  of  magnitude:  Re  =  1  x  103,  5  x  103,  5  x  104,  and 
5  x  105.  Two  aspect  ratios,  3.126  and  1.042,  were  employed,  pertinent  to 
apparatus  currently  available  for  experiments.  A  single  value  of  e  was  con¬ 
sidered,  namely  e  =  0.02,  which  lies  in  the  region  of  maximum  yaw  growth  rate 
occurring  in  BRL  gyroscope  experiments;  and,  finally,  t  =  0. 
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Our  interest  lies  primarily  in  the  side-plane  moment,  which  is  the  compo¬ 
nent  that  affects  the  yaw  of  the  projectile.  Figures  4  and  5  show  plots  of 
^LSM  versus  nutational  frequency,  x.  The  decrease  of  amplitude  and  flattening 

of  peak  with  decreasing  Reynolds  number  is  clearly  evident  for  both  aspect 
ratios.  Switching  the  aspect  ratio  from  3.126  to  1.042  changes  the  natural 
mode  primarily  being  excited  by  the  same  range  of  forcing  frequencies,  and  the 
new  amplitudes  consequently  differ  appreciably  from  the  old. 

For  comparison  of  our  output  with  other  results  we  look  at  computations 
by  the  Reference  6  method  (pressure  moments  only)  for  the  same  cases.  Corre¬ 
sponding  curves  lie  very  close  together  for  Re  =  5  x  lO4  and  5  x  105  for  both 
aspect  ratios,  and  divergence  increases  as  Reynolds  number  decreases. 
Generally,  the  Reference  6  method  gives  higher  peaks  than  the  present  method. 

Additional  comparisons  are  made  for  the  Clj^,  the  in-plane  moment  coeffi¬ 
cient,  in  Figures  6  and  7.  Again,  the  agreement  varies  in  an  inverse  manner 
with  Reynolds  number. 


VI.  YAW  GROWTH  RATE 

If  the  forcing  moment  of  Eq.  (2)  has  the  form 


MLY  +  i  mLZ  =  const  x  *>» 


then  the  following  expression  is  a  solution  to  Eq.  (2): 


i  =  K0eex^t  eiT^  -  K0eif^. 


(61) 


This  is  exactly  the  motion  assumed  in  Eq.  (3)  for  the  purpose  of  calculating 
the  liquid  moment.  The  amplitude  variation  is  governed  by  the  yaw  growth 
rate,  ex|.  The  growth  rate  factor,  e,  gives  a  measure  of  the  yaw  growth  or 

decay  in  one  nutational  cycle.  If  e>0,  |{j|  grows  until  the  small  yaw  assump¬ 
tion  is  no  longer  valid. 

Eq.  (61)  is  a  solution  to  Eq.  (2)  only  for  a  restricted  set  of  f's;  these 
values  are  the  solutions  to  the  functional  equation  obtained  by  substituting 
Eq.  (61)  into  Eq.  (2) : 


Iy  f2  -  Ix  f  -  Iy  M/I2  =  Pa5  F(f )  ,  (62) 

where  F  is  also  a  function  of  the  properties  of  the  fluid,  the  dimensions  of 
the  cavity,  and  the  spin-rate. 
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According  to  Stewartson-Wedemeyer,  F  of  Eq.  (62)  is  negligibly  small 
except  near  resonance.  A  resonance  condition  will  generally  occur  when 

Tn  ~  CR,  where  Tn  $  is  the  nutational  frequency  of  the  empty  shell 

and  Cr  is  a  natural  inertial  frequency  of  the  rotating  liquid. 

tn  =  £lx  +  Ox2  +  4  ly  M/$2) 

For  t  ~  CR,  F  can  be  approximated 
a  function  with  a  simple  pole1’3 

F  a  D  /  (f  -  Cr), 


1/2 


(2  Iv). 


by  the  first  term  of  the  Laurent  series  of 


where  the  residue,  D,  depends  on  the  parameters  of  the  problem.  In  this  case 
Eq.  (62)  yields  three  complex  solutions  for  f.  Only  those  solutions  are 

applicable  for  which  t  ~  CR;  generally  one  of  these  will  have  a  positive  e. 
Experiments  have  consistently  shown  that  for  a  significant  span  of  time,  the 
motion  is  described  by  an  exponential  solution  with  positive  e. 

In  our  analysis,  applicable  also  away  from  resonance,  the  K  of  Eq.  (61) 
and  the  moment  of  Eq.  (5)  are  substituted  into  Eq.  (2)  to  yield  the  following 
functional  equation  for  f: 


Iy  f2  -  Ix  f  -  Iy  M/$2  =  -i  (2tt  pa4c)  x  CLM  (f;  Re,  c/a). 


(63) 


This  equation  must  be  solved  by  iteration;  outputs  from  the  Stewartson- 
Wedemeyer  theory  furnish  guidelines  for  choosing  initial  estimates  of  f. 

The  motion  parameter  f  (=  x  -iex),  has  been  measured  in  gyroscope  experi¬ 
ments*.  The  calculations  of  Murphy6  indicate  that  viscous  shear  contributes 
significantly  to  the  liquid  moment  and  should  not  be  omitted.  Since  our 
pressure  moment  calculations  constitute  results  of  incomplete  theory,  we 


Applicability  of  gyroscope  experiments  to  simulation  of  projectile  angular 
motion  is  discussed  in  Section  2-7  of  Reference  7 . 
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shall  not  make  any  systematic  comparison  with  data  from  gyroscope  experi¬ 
ments.  However,  we  do  demonstrate  in  Figure  8  the  capability  of  our  program 
to  produce  yaw  growth  rate.  We  choose  two  cases;  the  first  was  treated  in 
Figure  8  of  Reference  4,  and  the  second  is  one  of  six  based  on  measurements 
described  in  Reference  10.  The  six  above-mentioned  cases  are  low  Reynolds 
number  experiments,  and  the  only  one  with  Re>104  is  shown  in  Figure  8. 

a 

The  parameter  that  is  varied  is  Tn  =  lx/ly  (since  M  =  0);  in  the  first 

cases  Iy  is  held  constant,  and  in  the  second  case  Ix  is  held  constant.  The 

quantities  t  and  e  are  measured  simultaneously;  data  are  presented  in  the  form 

of  gt  vs  t  plots.  Values  of  tn  are  chosen  that  lie  close  to  first  radial  mode 
eigenfrequencies  of  the  liquid,  namely,  k  =  3  for  c/a  =  3.15,  and  k  =  1  for 
c/a  =  1.04. 

Also  shown  in  Figure  8  are  yaw  growth  rates  predicted  by  the  Stewartson- 
Wedemeyer  theory,  which  likewise  treats  only  the  pressure  moment.  The  complex 
f  is  computed  from  the  cubic  equation  described  in  the  discussion  following 
Eq.  (63).  The  pairs  of  theoretical  curves  agree  very  well  for  Re  =  5.2  x  105 
but  differ  significantly  for  the  Re  =  1.24  x  104  case.  For  both  theories  the 
peaks  lie  to  the  left  of  the  experimental  peaks. 


VII.  DISCUSSION 

We  have  presented  a  detailed  treatment  of  the  linearized  problem  of  flow 
of  a  liquid  in  a  filled  spinning  cylinder  executing  angular  motion  at  small 
yaw,  and  we  have  exhibited  the  consequent  moments  exerted  by  the  liquid  on  the 
cylinder  walls.  The  angular  motion  takes  the  form  of  nutation  at  constant 
frequency  about  an  axial  point,  and  yaw  growth  at  an  exponential  rate;  the 
liquid  is  in  solid-body  rotation  before  the  angular  motion  begins.  This 
motion  approximates  the  actual  motion  of  the  liquid-filled  projectile  over  a 
significant  span  of  its  flight  history. 

The  emphasis  of  this  report  is  directed  primarily  to  the  method  of  deter¬ 
mining  output  rather  than  to  the  output  itself,  which  can  be  treated  in  more 
detail  in  later  studies.  We  feel  it  worthwhile  to  include  details  of  equa¬ 
tions,  formulas,  assumptions,  and  derivations  involved  in  the  lengthy  computa¬ 
tions  required  to  attain  the  liquid  pressure  and  moment  results. 

The  agreement  of  our  results  with  those  from  the  concurrent  work  of 
Murphy6,  who  extended  the  Stewartson-Wedemeyer  theory,  tends  to  support  the 
conclusions  from  both;  the  basic  assumptions  of  linearity  and  solid  body  are 
common  to  both  approaches.  In  addition,  both  make  use  of  boundary  condition 
corrections  which  are  applicable  for  large  Reynolds  numbers. 

Future  theoretical  work  will  include  the  additional  calculation  of  vis¬ 
cous  shear  moments.  Experimental  work  at  high  Reynolds  numbers  will  also  have 
to  be  performed  in  order  to  provide  an  adequate  data  base  for  validating  the 
theory. 


10.  W.  P.  D'Amico ,  Jr.,  and  T.  H.  Rogers,  "Yaw  Instabilities  Produced  by  Rapidly 
Rotatingx  Highly  Viscous  Liquids. ,"  AIAA  Paper  AIAA-81-0224,  AIAA  19th  Aero¬ 
space  Spences  Meeting ,  St.  Louis ,  Missouri,  January  12-15 ,  1981. 
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Figure  1.  Diagrams  of  Coordinates  and  Cylinder 


(b) 


Figure  2.  Comparison  of  Theoretical  and  Experimental  Endwall  Pressure 
Coefficients  (r=0.67,  c/a=3.148,  e=0.0). 
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Figure  3.  Comparison  of  Theoretical  and  Experimental  Endwall  Press 
Coefficients  (r=0.67,  c/a=1.051,  e=0.0). 


Figure  4.  Side  Moment  Coefficient:  Comparison  of  Results  of  Present 
Method  and  Method  of  Reference  6  (c/a=3.126,  e=0.02). 
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(c) 


Figure  4.  Side  Moment  Coefficient:  Comparison  of  Results  of  Present 
Method  and  Method  of  Reference  6(c/a=3.126,  e=0.02  (cont.) 


Figure  6.  In-Plane  Moment  Coefficient:  Comparison  of  Results  of  Present 
Method  and  Method  of  Reference  6  (c/a=3.126,  e=0.02). 


Figure  7.  In-Plane  Moment  Coefficient:  Comparison  of  Results  of  Present 
Method  and  Method  of  Reference  6  (c/a=1.042,  e=0.02). 
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LIST  OF  SYMBOLS 


ak 


c 


cross-sectional  radius  of  cylinder  [cm] 

least-squares  coefficients  for  series  of  axial 
eigenfunctions,  Eq.  (42) 

biorthogonal  coefficients  for  series  of  axial 
eigenfunctions,  Eq.  (40) 

half-height  of  cylinder 
=  c/a,  aspect  ratio 

perturbation  pressure  coefficient,  Eq.  (51) 
natural  oscillation  decay-rate  of  rotating  liquid/<j> 


natural  oscillation  frequency  of  rotating  liquid/^ 


CLIM 

CLSM 

Dp 

f 


I 

y 

k 


Ki 


liquid  moment  coefficient  =  CLSM  +  i  C^,  Eq.  (5) 
liquid  in-plane  moment  coefficient,  Eqs.  (6)  and  (60) 

liquid  side  moment  coefficient,  Eqs.  (6)  and  (60) 

disturbance  pressure  =  p-(£)>2,  Eq.  (48),  [pressure/ 
pa24>2] 

=(l-ie)x,  complex  representation  of  angular  motion 
Eq.  (4) 

moment  of  inertia  of  empty  shell  about  its  longi¬ 
tudinal  axis 

transverse  moment  of  inertia  of  empty  shell  about 
its  center  of  gravity 

index  of  axial  eigenfunction  and  eigenvalue,  Eqs. 

(28)  and  (31) 

yaw  amplitude  at  time  t  =  0 

_  KoeET<flt,  yaw  amplitude  at  time  t,  Eqs.  (3)  and  (4) 
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i 

x  (and  x)  coordinate  of  pivot  point 

mL 

mass  of  liquid  in  cylinder  =  2irpa2c  [g] 

MLY’  mlz 

y  and  z  components,  respectively  of  liquid  moment 
[g  cm2/s2] 

MLZS’  mlze 

sidewall  and  endwall  contributions,  respectively, 
to 

MLZT’  mlzb 

top  and  bottom  wall  contributions,  respectively,  to 

\z 

Mp  M2 

non-dimensional  liquid  moment  coefficients,  in-plane 
and  side  plane,  respectively,  Eq.  (58) 

M1E*  M2E 

endwall  contributions  to  M,  and  M9,  Eqs.  (56),  (57c), 
(57d),  (59) 

Mls*  M2s 

sidewall  contributions  to  M,  and  M9,  Eqs.  (56),  (57a), 
(57b),  (59) 

A 

M 

aerodynamic  (or  gravity)  moment  parameter,  Eq.  (2) 

fs‘2] 

nYE’  nZE 

components  of  the  projection  in  the  y,  z  plane  of  a 
unit  vector  along  projectile  axis 

P 

pressure/(pa24>2) 

p  =  pr  +  1  Pi 

non-dimensional  r,  x  variation  of  perturbation 

pressure,  Eq.  (17d),  [pressure/(K  pa2^2)] 

0 

B  H 

series  solution  contribution  to  p,  Eq.  (28) 

-pp 

particular  solution  contribution  to  p,  Eq.  (22) 

★ 

p 

perturbation  pressure/(Kopa2<t>2) ,  Eq.  (10) 

<  CL 

coefficient  of  sin  X^x  in  p^  series,  Eq.  (28) 

P 

unperturbed  pressure/(p a 2 4>2 ) ,  Eqs.  (10)  and  (11) 

Pi.  P2 

perturbation  pressure  coefficients,  Eq.  (50), 
[pressure/ (KQpa 2  412 )] 

r 

radial  coordinate  in  inertial  system/a 
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IE  IE 


r 

Ri »  •  •  •  *  R** 


Re 

t 

t 


radial  coordinate  in  non-rotating  aeroballistic 
system/a 

functions  of  radial  variation  in  separated  flow 
solutions  uH>  yH>  wH,  pH,  respectively,  Eqs.  (26) 

and  (C.l ) 

Reynolds  number  =  a2  4»/ v 
time  [s] 

=  <j>t 


u,  v,  w  non-dimensional  radial,  azimuthal,  and  axial  velocity 

components,  respectively,  in  inertial  system  [velo¬ 
city/^)] 

u,  v,  w  non-dimensional  x,  r  variation  of  perturbation  velo¬ 

city  components,  Eq.  (17)  [velocity/(Ko<fia)] 

H’  -H’  -H  series  solution  contribution  to  u,  v,  w,  Eq.  (28) 


*  *  * 
U,  V,  w 

U,  V,  w 


★  ★  ★ 


Uf 


V, 


w 


u 


k’ 


U,  V,  W 


X 


i  » 


•  •  •  J 


X4 


x,  y,  z 

x,  y>  z 


particular  solution  contribution  to  u,  v,  w,  Eq.  (22) 

non-dimensional  perturbation  velocity  components  in 
inertial  system,  Eqs.  (10)  and  (Al)  [velocity/(Ko<fa)] 

non-dimensional  radial,  azimuthal,  and  axial  velocity 
components ,  respectively ,  in  aeroballistic  system 
[velocity/(<£a)] 

non-dimensional  perturbation  velocity  components  in 
aeroballistic  system,  Eq.  (A.l)  [velocity/U^a)] 

coefficients  of  sinx^x  and  cosx^x  in  u^,  v^, 
series,  Eq.  (28) 

non-dimensional  radial,  azimuthal,  and  axial  velocity 
components  of  unperturbed  flow,  Eq.  (11)  [velocity/ 

(4>a)l 

functions  of  axial  variation  in  separated  flow  solu¬ 
tions  U|_j ,  yH,  Wj_| ,  pH,  respectively,  Eqs.  (26)  and  (Cl) 

non-dimensional  rectangular  coordinates  in  inertial 
system  (x-axis  along  trajectory)  [length/a] 

non-dimensional  rectangular  coordinates  in  aeroballist¬ 
ic  system  (x-axis  along  cylinder  axis)  [length/a] 
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y  i»  •  •  •  y  6 


functions  describing  radial  variation  of  perturbation 
flow  variables,  Eqs.  (34)  and  (35) 


e 


6c 


E 


0, 


e 


x 

xl 

V 

'"W 

P 

T 

f 


angle  in  vertical  plane  measured  from  the  x-axis  to 
the  velocity  vector 

angle  in  horizontal  plane  measured  from  the  x-axis  to 
the  velocity  vector 

correction  term  in  endwall  boundary  condition,  Eqs. 
(29)  and  (30) 

=  (1/t)  *  yaw  growth  per  radian  of  nutation 

polar  angles  in  inertial  and  aeroballistic  systems, 
respectively,  Eqs.  (7)  and  (8) 

separation  of  variables  constant,  Eq.  (27) 

eigenvalue  in  the  axial  problem,  Eqs.  (28)  and  (31) 

kinematic  viscosity  of  liquid  [cm2/s] 

vector  describing  angular  motion  of  cylinder,  Eq.  (1) 

density  of  liquid  [g/cm3] 

nutational  frequency  of  cylinder/^ 

=  T<f>t,  angular  orientation  of  x-axis  in  the  x,  y,  z 
system 

spin  rate  of  cylinder  [rad/s],  taken  to  be  positive 
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APPENDIX  A.  DERIVATION  OF  WALL  BOUNDARY  CONDITIONS 


The  wall  boundary 
*  *  * 

u,  v,  w  at  r  =  1  and  x 


(11), 


conditions,  Eq.  (15),  must  be  restated  in  terms  of 
=  ±c .  Recalling  the  definitions  of  Eqs.  (10)  and 


_  ★ 
u  =  dr/dt  =  -K  u, 

O 


v  =  r  de/dt  = 


r-K  v, 

o 


w  =  dx/dt  =  -K  w 

O 


(Al) 


u  =  dr/dt 


■K  u, 

O 


v  =  r  de/dt  = 


r-K  v, 

O 


w  =  dx/dt  =  -K  w, 

O 


where  t  e  4>t. 

The  velocity  transformation  between  earth-fixed  and  aeroball istic  systems 
is  obtained  by  differentiating  Eq.  (9)  with  respect  to  time  and  substitut¬ 
ing  the  expressions  of  Eq.  (A  1)  for  the  derivatives: 

★  • 

u  =  u-(x-£)eeT<t>t[eT  cos (<j>  -e)  +  ( 1  -t )  sin(<j>  -e)]  +  0(K  ) 

1  1  o 

•  * 

v  =  v-(x-«,)e£T4>t[-(l-T)  cos(<j>  -e)  +  ex  sinU  -e)]  +  0(K  )  /AO, 

1  1  o  \HC.) 

w  =  w-r  eET<’>t[-ex  cos(<j)  -e)  -  (1-x)  sin(<j>  -e)]  +  0(K  )  . 

1  1  o 


The  tilde  superscripts  can  be  dropped  from  the  bracketed  terms  of  Eq.  (A2) 
without  changing  the  order  of  error. 

★ 

We  now  evaluate  Eq.  (A2)  at  r  =  1 ,  applying  the  conditions  u(r=l)  = 

*  * 

v(r=l )  =  w(r=l )  =  0.  Thus, 


u(r=l ) 


u(r=l)  +  (rr=1-l)  (9u/9r)~=1 

•  * 

0  +  K  e£Tl^t(x-£)  cosU  -e)  (9u/ar)~  , 

O  1  I  —  l 


(A3) 


*  * 

Similar  expressions  apply  to  v(r=l)  and  w(r=l).  The  coefficients  of„the 
first  order  terms  (in  K  )  contain  gradients  which  might  be  large  at  r  =  1 

O 

when  the  Reynolds  number  is  large.  However,  at  this  stage  we  are  only 
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considering  perturbations  from  the  point  of  view  of  the  small  parameter  K  . 
Hence,  to  zeroth  order,  0 

*  *  * 

u(r=l)  =  v(r=l)  =  w(r=l)  =  0,  and  consequently 
u(r=l )  =  -(x-£)eeT<*>t[  ex  cos(4>  -e)  +  (1-x)  s i n ( -e)] 

(A4) 

*(r=l )  =  -(x^)eET^[-(l-T)  cos ( ij>  -e)  +  ex  sin(<j>  -e)] 

i  i 

v(r=l)  =  -eGT<^)t[-ex  cos(ij)  -e)  -  (1-x)  sin($  -e)]. 

i  i 

Similar  consideration  applied  to  the  endwalls  leads  to  the  conclusion 
?-*_*_ 

that  u(x=±c)  =  v(x=±c)  =  w(x=±c)  =  0  +  0(K  );  Eqs.  (A2)  now  yield 

o 

u(x=±c)  =  -(±c-£)eET<I>t[+ex  cos ( 4>  -0)  +  (1-x)  sin(4>  -e)] 

i  i 

v(x=±c)  =  -(±c-2.)eETI^>t[-(l-x)  cos($  -e)  +  ex  sin(<j>  -e)]  (A5) 

i  i 

w(x=±c)  =  -r  eex<f>t[-ex  cos(<j>  -e)  -  (1-x)  sinU  -e)]. 

i  i 
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APPENDIX  B.  AXIAL  BOUNDARY  CONDITIONS 

The  first  of  Eqs.  (20),  namely,  u-iy  =  0  at  r=0,  is  obtained  merely  by 
setting  r=0  in  Eq.  (18a).  The  next  condition,  namely,  w(r=0)  =  0,  is  ob¬ 
tained  by  multiplying  Eq.  (18d)  by  rz,  then  setting  r=0. 

If  one  differentiates  Eq.  (18a)  with  respect  to  r,  then  sets  r  =  0 
(remembering  that  wx(r=0)  =  0  because  w(r=0)  =  0),  one  obtains 

2ur(r=0)-ivr(r=0)  =  0.  (Bl) 

If  one  multiplies  Eq.  (18c)  by  r2,  differentiates  with  respect  to  r, 
then  sets  r=0,  one  obtains 

ip(r=0)+(l/Re)[-yr(r=0)-2iur(r=0)]  =  0.  (B2) 

Application  of  Eq.  (Bl)  here  leads  to  the  third  condition,  namely, 

p(r=0)  =  0.  (B3) 


Alternatively,  the  boundary  conditions  at  r=0  can  be  derived  only  on 
the  basis  of  continuity  and  single-valuedness. 5 
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APPENDIX  C.  SEPARATION  OF  VARIABLES  IN  THE 
LINEARIZED  NAVIER-STOKES  EQUATIONS 

Our  perturbation  equations  follow  directly  from  the  linearized  Navier- 
Stokes  equations.  After  the  t  and  e  variations  are  removed,  as  in  Eqs.  (17), 
the  perturbation  equations  are  in  the  form  of  Eqs.  (18).  The  form  of  the 
solution  to  these  equations  was  obtained  by  separation  of  variables.  Since 
this  procedure  cannot  be  carried  out  strictly  in  the  usual  manner,  it  is  out¬ 
lined  here. 

The  variables  are  separated  as  in  Eq.  (26)  and  substituted  into  Eqs. 
(18).  The  manipulations  will  be  outlined  only  for  Eq.  (18b)  since  the  others 
are  treated  in  the  same  way.  After  collecting  terms,  Eq.  (18b)  becomes 


[Re-  (Rn  +  R'/r  -  2R  /r2 )  -  i(f-l)]  X  + 

ill  i 

Re-1R  X"  +  2  [Re-1  i/r2  +  1]  R  X  -  R'X  =  0, 

1  1  2  2  4  4 


(Cl) 


where  the  prime  indicates  differentiation  with  respect  to  the  argument.  It 
is  convenient  to  denote  each  coefficient  of  the  functions  of  x  by  a  single 
symbol.  Thus  Eq,  (Cl)  is  rewritten  as 

g  (r)  X  +  g  (r)  X"  +  g  (r)  X  -  g  (r)  X  =0.  (C2) 

1  1  2  13  21*  4 


Divide  this  by  g  and  differentiate  with  respect  to  r  to  obtain 

4 


(g  /9  )'  x  +  (g  /g  )'  X"  +  (g  /g  )'  X  =0. 

141  2  4  1  34  2 

Next,  divide  by  (g  /g  )'  and  differentiate  with  respect  to  r.  This  gives 

3  4 

[(g  /g  ) 7(g  /g  )']'  x  +  [(g  /g  ) 7(g  /g  )']'  X"  =  0.  (C3) 

14  34  1  24  31*  1 

The  usual  argument  employed  in  separating  variables  shows  that 

X"  =  const.  X  , 

1  1 


which  is  the  form  given  in  Eq.  (27). 

The  same  process  applied  to  Eqs.  (18d)  and  (18a)  yields 


X"  =  const.  X 
3  3 
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and 


X  =  const.  X  , 

2  i 

respectively.  When  the  above  three  results  are  substituted  into  Eq.  (C2), 
the  result  is 


X  =  const.  X  . 

b  1 

The  boundary  conditions  dictate  the  solution  to  the  harmonic  equations,  as 
shown  in  Eq.  (28).  The  particular  solution  stated  in  Eqs.  (22)  is  obtained 
by  setting  equal  to  zero  the  constants  of  X"/X  =  constant. 

i  i 
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APPENDIX  D.  THE  ENDWALL  BOUNDARY  CONDITION 


After  various  transformations,  the  flow  problem  to  be  solved  is  given  by 
Eq.  (18),  the  linearized  Navier- Stokes  equation,  plus  the  boundary  conditions 
on  the  axis,  Eq.  (25),  on  the  sidewall,  Eq.  (23),  and  the  endwalls,  Eq.  (24). 
The  method  of  solution  adopted  is  separation  of  variables  in  order  to  get  the 
normal  modes.  The  axial  and  sidewall  boundary  conditions  can  be  satisfied 
with  the  modal  solution,  but  only  the  condition  of  no  flow  through  the  end 
wall,  Eq.(24c),  can  be  satisfied  there.  The  remaining  no  slip  conditions,  Eq. 
(24a,  b),  cannot  be  satisfied.  Rather  than  the  viscous  boundary  conditions, 
we  satisfy  the  condition  w^  =  0  that  would  be  imposed  for  an  inviscid  fluid 

and  accept  the  values  of  u^  and  v^  on  the  endwalls  from  the  solution,  as  in  an 

inviscid  flow. 

The  proper  normal  modes  are  obtained  only  for  \  real  in  Eq.  (27).  How¬ 
ever,  the  modal  solution  is  inconsistent  with  the  linearized  Navier-Stokes 
equations  and  the  no  slip  boundary  conditions.  Since  we  have  had  to  drop  two 
of  the  three  endwall  conditions,  we  can  conjecture  that  a  boundary  layer  must 
be  inserted  in  order  to  satisfy  these.  This  is  what  Wedemeyer3  did  to  correct 
the  inviscid  solution  of  Stewartson1.  This  is  the  type  of  correction  needed 
here.  It  is  done  in  the  same  spirit  as  Wedemeyer1 s  correction  but  must  be  de¬ 
rived  as  a  correction  to  Eqs.  (18)  which  are  not  Stewartson' s  inviscid  equa¬ 
tions. 

The  technique  of  matched  asymptotic  expansions  is  used  here  since  it  pro¬ 
vides  a  systematic  method  of  handling  this  type  of  singular  perturbation  prob¬ 
lem;  the  results  are  more  easily  interpreted  than  in  Wedemeyer's  approach,  and 
some  difficulties  in  his  approach  are  avoided.  The  technique  is  described  by 
Van  Dyke11.  However,  the  way  the  technique  is  applied  is  tailored  to  the 
particular  problem  we  solve. 


To  avoid  a  proliferation  of  indices,  the  notation  is  changed  for  this 


Appendix 

only.  The 

sub-bar  and  sub-H 

are 

dropped. 

see  Eq.  (21),  and  we 

let  C  = 

c.  The  flow 

variables  are  now  u. 

v,  w,  and 

p.  The  linear  Eqs. 

(18) 

are  written  with  the 

linear  operator. 

L(u ,  v,  w. 

P)  : 

=  0- 

(Dl) 

Only  the 

boundary  conditions  at  x  =  C 

are 

needed  in 

this  derivation: 

u  =  ig 

v  =  g 

at  x  =  C 

(D2) 

w  =  0 

g  =  -  [2f  (i-f)/(i+f)]  c 

Van  Dyke,  Perturbation  Methods  in  Fluid  Mechanics ,  Academic  Press, 
New  York ,  N.YTJ~TUW. - “ - 
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The  asymptotic  solution  to  Eq.  (Dl)  with  Eq.  (D2)  for  Re  +  »  is  obtained 
by  constructing  inner  and  outer  solutions  with  appropriately  scaled  inner  and 
outer  variables.  The  outer  variables  are  unchanged,  which  conforms  with  Eq. 
(18).  The  outer  solution  has  the  asymptotic  form 


u  ~  u0 

(x. 

r. 

Re) 

+  A1 

(Re) 

uol 

(x. 

r. 

Re)  +  ... 

< 

l 

< 

0 

(x. 

r. 

Re) 

+  a2 

(Re) 

vol 

(x. 

r. 

Re)  +  ... 

w  ~  w 

0 

(x. 

r. 

Re) 

+  A3 

(Re) 

wol 

(x. 

r. 

Re)  +  ... 

p  ~  p0 

(x. 

r. 

Re) 

+  A4 

(Re) 

Pol 

(x. 

r. 

Re)  +  ... 

where  A.  -*■ 

0  as 

Re 

“1 

,  and 

this 

form 

is 

valid  for  Re 

fixed.  For  the  inner  solution  x  is  replaced  by  y  =  C  -  x. 


(D3) 


with  x  and  r 

and  inner  variables 


y  =  y/e  (Re)  (D4) 

r  =  r 

are  introduced  where  e  +  0  as  Re  >  «;  in  this  Appendix  super-bar  does  not 
denote  complex  conjugate.  The  inner  solution  has  the  asymptotic  form 


u  ~  u.  (y,  F)  +  6n  un  (y,  r)  +  ... 

v  ~  v.  (y,  r)  +  «21  vn  (y,  r)  +  ... 

w  ~  63  w.  (y,  r)  +  631  wn  (y,  r)  +  ... 

P  ~  «4  Pi  (y.  r)  +  641  pn  (y,  r)  +  ... 


(D5) 


where  6.  +  0  as  Re  +  <5.,  <<  6.,  and  the  form  is  valid  for  Re  +  00  with  y 

J  J  ^  J 

and  r  fixed. 

The  scales  A.,  6.,  6.,  and  e  must  be  determined  from  matching  the  inner 
J  J  Jl 

and  outer  solutions  and  from  whatever  information  the  problem  provides.  How¬ 
ever,  the  terms  in  Eq.  (03),  uQ  (x,  r,  Re),  etc.,  and  uQl  (x,  r,  Re),  etc., 

are  not  in  proper  form  for  matching  because  they  are  functions  of  Re.  This 
follows  from  the  fact  that  the  first  terms  satisfy 
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(D6) 


and  similarly  for  the  second  terms  u^,  etc.  Matching  requires  that  the  only 
dependence  on  Re  be  in  the  scales  A..  The  terms  of  Eq.  (D3)  can  also  be  ex- 

•J 

panded  in  asymptotic  series: 


uo  ~  uo°  r)  +  Re_1  uq1  ^x*  + 


(D7) 


etc. 


and 


uoi~  uoi°  (x*  r‘)  +  Re_1  u02 1  (x»  r)  +  •••*  (D8) 

the  form  of  the  expansion  being  determined  by  Eq.  (D6).  The  functions  uQ° 
etc.,  and  u  0  ,  etc,  satisfy  the  inviscid  equations  and  are  not  functions  of 

1  U1  i 

Re;  ,  etc.,  and  uQ^  ,  etc.,  satisfy  a  non-homogeneous  form  of  the  inviscid 

equations.  This  explains  why  we  can  get  a  solution  with  only  the  one  boundary 
condition  on  normal  velocity.  The  process  described  in  this  paragraph  is 
tailored  to  the  problem  we  solve;  a  more  direct  approach  would  normally  be 
used. 


The  asymptotic  matching  principle11  can  now  be  applied.  Only  the  results 
are  given: 


e  =  Re 


-1/2 


51  =  ^2  "  64  "  1 


A.  =  Re"1/2,  j  =  1,  2,  3,  4 

•J 


53  "  Re 


-1/2 


The  first  terms  in  Eq.  (D5)  satisfy  the  boundary  layer  equations: 

i  (f-l)u.  -2v.  =-p.  -+u.  -  - 

i  i  r  i  y  y 

i  (f  -  1)  v.  +  2Ui  =  i  p./F  +  v.  -  - 


(D9) 


p .  -  =  0 

y 


(r  Uj)  p  -  i  v.  -  F  w.  -  =  0 
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with  boundary  conditions  on  the  endwall ,  obtained  from  Eq.  (D2), 


ui  (0,  r)  =  ig 

v.  (0,  F)  =  g 

wi  (0,  r)  =  0 

and  at  the  edge  of  the  boundary  layer,  obtained  from  matching, 

Ui  («,  r)  =  uo°  (C,  r) 

vi  (°°»  ?)  =  vo°  (C*  r) 

Pi  (-.  ?)  =  p0°(C,  r) 

or  more  generally  p^(y,  r)  =  p°  (C,  r).  No  condition  on  («>,  r)  is  obtained 

or  allowed.  The  solution  to  Eq.  (D9)  is  easily  obtained  and  is  required  in 
the  matching. 

The  boundary  condition,  Eq.  (29a),  must  be  interpreted  as  a  condition  on 
the  outer  solution  which  is  given  by 

w  ~wo°  (x,  r)  +  Re'1/2  wol°  (x,  r)  +  0  (Re-1)  (DIO) 

The  functions  here  satisfy  only  one  boundary  condition  on  the  wall: 


wQ°  (C,  r)  =  0 


wol°  r)  =  Rel/2  6c 


x  =  C 


(Dll) 


the  latter  being  obtained  by  matching.  Using  Eqs.  (Dll)  and  (DIO)  we  obtain 


w  -  <5cwx  =  0  +  0  (Re"^)  at  x  =  C, 


( D12) 
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which  is  Eq.  (29a).  This  development  shows  that,  rather  than  solving  Eq. 

(01),  we  could  solve  two  inviscid  type  problems  with  boundary  conditions  Eq. 
(Dll)  and  combine  them  linearly  to  get  Eq.  (DIO).  Actually  this  is  not  prac¬ 
tical  because  the  radial  variation  and  no  slip  conditions  at  the  sidewall  must 
be  found. 

The  above  asymptotic  solution  was  obtained  using  two  terms  in  the  outer 
solution  and  one  term  in  the  inner  solution.  To  obtain  the  next  approximation 
to  the  boundary  condition,  three  terms  in  the  outer  and  two  terms  in  the  inner 

-3/2 

would  be  necessary.  The  boundary  condition  will  have  an  error  0  (Re  )  and 
should  allow  a  more  accurate  solution  for  lower  Re. 
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APPENDIX  E.  THE  B I ORTHOGONAL  EXPANSION 

Separation  of  variables  leads  to  the  result  that  the  X  in  Eq.  (26) 

i 

satisfysthe  harmonic  equation,  Eq.  (27),  where  X  is  a  separation  constant. 
Considering  w  and  applying  the  boundary  conditions  (29)  yield  the  following 
H 

boundary  value  problem,  assuming  f,=0, 

Z"  +  \z  1  =  0 


Z'(-c)  -  x2sc  Z(-c)  =0  (El) 

V  (c)  +  X26c  Z(c)  =  0 

where  X  (x)  is  replaced  by  Z'(x)  because  we  wish  to  develope  the  expansion 

3 

for  x  needed  in  Eq.  (38).  The  eigenvalue,  X,  is  determined  from 

(1 -6c2X 2 )  sin  2X  c  =  2cX  cos  2x  c  (E2) 

J  j  j  j 


which  can  be  split  into  two  equations,  one  of  which  is  Eq.  (31)  and  is  the 
appropriate  one  for  our  problem.  A  denumerable  set  of  complex  solutions, 

X  ,  exists.  Even  though  the  differential  equation  in  (El)  is  self-adjoint, 
j 

the  system  (El)  is  not,  which  means  that  complex  eigenvalues  are  to  be  ex¬ 
pected.  The  non  self-adjoint  problem  and  the  presence  of  the  eigenvalue  in 
the  boundary  conditions  make  this  a  nonstandard  problem.  The  eigenfunctions 
are  not  orthogonal  and  several  of  the  usual  results  for  expansion  of  a 
function  in  a  series  of  eigenfunctions  do  not  apply. 

This  type  of  problem  can  be  attacked  by  introducing  the  adjoint  problem: 


Y"  +  XZY  =  0 


_  2 


Y'(-c)  -X  6c  Y(-c)  =  0 


Y'(c)  +  X26c  Y(c)  =  0 


(E3) 


where  -  denotes  complex  conjugate. 
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An  essential  part  of  the  theory  for  the  expansion  of  a  function  in  a 
series  of  eigenfunctions  is  the  proper  definition  of  an  inner  product.  In 
this  problem  solutions  Z  (eigenvalue  x  )  to  Eq.  (El)  and  Y  (eigenvalue  A  ) 

j  J  k  k 

to  Eq.  (E3)  are  needed  to  define  the  inner  product  (Z  ,  Y  ).  It  can  be 

j  k 

shown  that  the  inner  product  must  be  defined  as 


(Z  ,  Y  ) 
j  k 


Z  Y  dx  -  Sc[Z  (c)  Y  (c)  +  Z  (-c)  Y  (-c)]. 
j  k  j  k  j  k 


The  boundary  terms  in  this  expression  arise  from  the  presence  of  A  in  the 
boundary  conditions.  The  functions  Z  and  Y  are  biorthogonal  if 

j  k 

(Z  ,  Y  )  =  0. 
j  k 


Any  solution  Z  is  biorthogonal  to  any  solution  Y  if  j  +  k.  However,  Z 
j  k  j 

and  Z  or  Y  and  Y  are  not  biorthogonal. 
k  j  k 


A  more  or  less  arbitrary  function  can 
functions,  Z  ,  (or  equivalently  Y  ) 


be  expanded  in  terms  of  the  eigen- 


where  C  =  (f,  Y  )/(Z  ,  Y  ) 

k  k  k  k 


with  all  the  usual  desirable  properties.  In  particular  the  function  x  can 
be  expanded  as  in  Eq.  (39)  with  the  coefficients  b  given  in  Eq.  (40). 

k 
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